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Detection and correction of phase jumps in a phase sequence 

Field of the invention 

The invention relates to a conununication system comprising at least a 
transnodtter and a receiver intended to receive symbols coming from a PSK modulation, and 
comprising estimation means for estimating a frequency error relating to a symbol based on a 
sequence of s3mibol phases. The invention also relates to a receiver intended to be used in 
5 such a commuoication system. 

The invention also relates to a method of estimating a frequency error relating 
to a received symbol, coming from a PSK modulation, based on a sequence of symbol 
phases, and to a method of detecting and correcting phase jumps in an initial phase sequence 
of symbols coming from a PSK modulation. 
10 The invention finally relates to computer programs comprising instmctions for 

implementing such methods. 

Backgroimd of the invention 

Such an estimation algorithm of a frequency error relating to a received 

IS symbol based on a phase sequence is described, for example, in paragraph 4.1 (page 107) of 
the article «Feedforward Frequency Estimation for PSK : a Tutorial Review» by 
M. Morelli and U. Mengali, published in the joiimal «Eiiropean Transactions on 
Telecommunications, vol. 9, no. 2, March ~ April 1998». This algorithm is known by the 
name of Tretter algorithm, or least squares method. 

20 To obtain such a phase sequence, it is known that a phase estimation algorithm 

is used, which estimates the phase relating to a received symbol on the basis of decisions 
made on various received symbols. For example, the expectation maximization algorithm is 
used, which is described in paragraph 3.3 of the conference report of the «Intemational 
Conference on Communications, New Orleans, USA, 1-5 May, 1994, vol. 2, pp. 940 and 

25 945», entitled «Comparison between digital recovery techniques in the presence of 
frequency shift» by F. Daffara and J. Lamour. 
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The problem posed is the following: in a ESK modulation comprising 2" 
points, two adjacent points have a phase difference of . When the frequency error is such 
that, based on a certain symbol, an error is made in the decision, this error is translated by a 
phase jump of in the sequence of phase estimates obtained. The number of phase 

jumps that may be obtained depends on the number of symbols contained in the packet and 
the initial frequency difference. 

When the sequence of phases, which are used for applying the Tretter 
algorithm, includes one or various phase jimips, the frequency estimate obtained is 
inaccurate. The invention notably has for its object to provide a solution to this problem. 

Summar y of the invention 

For this purpose, a communication system according to the invention and as 
described in the opening paragraph is characterized in that said receiver comprises 
calculation means for calculating a phase sequence, called initial sequence, based on 
decisions made on symbols, and means for detecting and correcting phase jumps in this 
initial sequence, to supply a phase sequence, called final sequence, to said frequraicy error 
estimation means. 

Advantageously, said means for detecting and correcting phase jumps 

comprise: 

- modifying means for modifying said initial sequence, so as to produce a plurality of 
modified sequences which each compensate for a pheise jump configuration, 

- calculation means for calculating straight Ime equations which determine the initial 
sequence and the modified sequences, 

- calculation means for calculating for the initial sequence and the modified sequences a 
mean difference between the initial or modified phases and the phases produced by the 
corresponding straight line equation, said final sequence being formed by the sequence 
whose mean difference is niiiumal. 

Brief description of the drawings 

These and other aspects of the invention are apparent from and will be 
elucidated, by way of non-limitative example, with reference to the embodiment(s) described 
hereinafter. 
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In the drawings: 

- Fig. 1 is a diagram of an example of a conmiunication system according to 

the invention, 

5 - Fig. 2 is a flow chart describing the steps of a phase estimation method, 

- Fig. 3 is a flow chart diagram describing the operations used by the phase 
jump correction and detection means according to the invention, 

- Fig. 4 is a representation of a phase sequence includmg a phase jump, 

- Fig. 5 is a representation of a phase sequence including two phase jumps, 
10 - Fig. 6 is a representation in the form of curves of the results obtained with 

the invention. 



Description of preferred embodiments 

In Fig. 1 is represented an example of a communication system according to 

15 the invention. This communication system comprises interactive user terminals 1, which are 
transmitters within the meaning of the invention, and a head-end station 2, which is a receiver 
within the meaning of the invention. The head-end station 2 transmits signals in a first 
fiequency band Ku (12-14 GHz). These signals are relayed to interactive user terminals 1 by 
a satellite 3. The interactive user terminals transmit signals in a second jfrequency band Ka 

20 (20-30 GHz). These signals are relayed to tiiie head-end station 2 by the satellite 3. 

Each terminal 1 comprises a data source 10 and channel coding means 12. The 
channel coding means deliver packets of N symbols which contain preamble symbols and 
data symbols. These packets are then transmitted to filter means 13 and, finally, to 
modulation means 14, which use a local oscillator that has a firequency fc. 

25 The head-end station 2 comprises demodulation means 20 which use a local 

oscillator that has a frequency t+Afo/Ts (where Afo is a normalized fi-equency difference 
relative to the symbol fi-equency, and Ts is the duration of the symbols), and an initial phase 
©0. The head-end station 2 also comprises filter means 21 and sampling means 22, which 
sample the output signal of the filter 21, to deUver symbols called received symbols. The 

30 received symbols are transmitted to frequency recovery means 23 which estimate the 

normalized firequency difference Mq and correct the received symbols to compensate for the 

estimated difference A/o . The frequency recovery means 23 deliver frequency-corrected 
symbols. These fi^quency-corrected symbols are transmitted to phase recovery means 24. 
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The phase-corrected and frequency-corrected symbols are finally transmitted to channel 
decoding means 25 which deliver data 27. 

In the example described here, the phase recovery means 24 comprise phase 
and frequency estimation means 50 and phase correction means 51 . The phase estunation 
means 50 are described with reference to Fig. 2. They are formed by a loop mtended to be 
passed through L times. In the foUowmg of the description the index m (m=l to L) is a loop 
counter. Each value of m thus corresponds to one loop path. The loop comprises: 

- conventional phase estimation means 52 (for example, expectation maximization) for 
producing an initial phase sequence SI relating to a symbol packet (m=l to L, and 

fc==l to q, where q<N), 

- means 53 for calculating a frequency error relating to said packet, based on the initial 
sequence SI, 

- means 54 for correctmg the frequency of the symbols of the packet, to correct said 
frequency error, 

- loop means 55 which supply the frequency-corrected symbols rk^"^^ to phase estimation 
means 52 for a following patii through the loop, 

- and loop output means 56 which supply to the phase correction means 51 the phases 
estimated by the phase estimation means 52 (or, directly, the correction to be made), and 
the symbols to be corrected. 

The calculation means 53 calculate a frequency error estimate Kfm relating to 
the symbols rA^'"*"^^ of a same packet, based on the initial phase sequence SI produced by the 
phase estimation means 52. The values of this phase sequence SI are between -oo and +oo. 
The calculation means 53 comprise: 

- means 60 for detecting and correcting phase jmnps to correct this initial sequence SI and to 
deliver a final sequence S2, 

- means 62 for calculating the slope of a straight line which is closest possible to the values of 
the sequence S2, by applying the Tretter algorithm. The slope obtained forms a frequency 
error estunate &fm relating to the symbols of the packet It is tiiis frequency error that is 
transmitted to the frequency correction means 54. The symbols obtained after firequency 
correction, rk^"^^ = rjS^'^^^ .e"*^^^ ^ , are transmitted to the phase estimation means 52 for a 
new path tiirough the loop. During the last path tiirough the loop (m=L), the symbols /-jj^"^^ 
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to be corrected and th6 phase* correctioii e ^ to be made to these symbols are 
transmitted to the phase correction means 51. 

The means 60 for detecting and correcting phase jimips are represented in Fig. 
3. They comprise: 

5 - means 100 for modifying the initial sequence to compensate for a plurality of phase jump 
configurations; the correction means 100 produce a plurality of modified sequences which 
correspond each to the correction of a phase jimip configuration; 

- means 110 for calculating straight line equations which determine the initial sequence and 
the modified sequences, 

10 - calculation means 120 for calculating for the initial sequence and the modified sequences a 
mean difference between the initial or modified phases and the phases produced by the 
corresponding straight line equation, said final sequence being formed by the sequence 
v^ose mean difference is minimal. 

The use of the means 60 for detecting and correcting phase jumps is different 

1 S and depends on the number of phase jxmxps one wishes to correct. But the method applied 

remains the same. The invention can thus be applied to any number of phase jumps. Now two 
examples of embodiment will be described of the means 60 for detecting and correcting 
phase jmnps: 

- a first example, in which the means 60 for detecting and correcting phase jumps are 
20 intended to correct a single phase jimip of ±^ ; 

* a second example, in which the means 60 for detecting and correcting phase jumps are 
intended to correct two phase jumps of ± y ia the same direction; this second example 

corresponds to the most probable case where the transmitted packets are ATM cells of S3 
octets. 

25 In the examples that will be described, the modification means 100 modify the 

initial sequence SI phase by phase. But, for diminishing the number of calculations to be 
performed, it is possible to modify the sequence phase-groi^) by phase-group. This is 
equivalent to compensating only for certain configurations of phase jumps. 

The first example of embodiment of the means 60 for detecting and correcting 

30 phase jumps is described with reference to Fig. 4. In Fig. 4 is represented an example of an 

initial sequence S 1 comprising a phase jump of + y . The means 60 for detecting and 
correcting phase jumps have for their Amotion to detect the position and the direction of the 
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phase jump and then correct it. Therefore, as indicated in Fig. 3, they perform the following 
operations: 

a) The Tretter algorithm is applied to the initial sequence SI formed by phases cpj Q=0 to q-1) 
to obtain the straight line equation Do which determines this sequmce. This equation is 
5 written as: 

Do: y=ao.x+bo with: 

=aS'->8.S and ^yS-fiS 



10 These expressions can easily be derived from the calculations shown on pages 523 and 524 
of the title «Numerical Recipes in C, the art of scientific computing, second edition» by 
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, published by Cambridge 
University Press in 1995, while considering that the uncertainty of the phases is constant 
whatever j. 

15 b) An initial mean difference (o-^ )^ is calculated between the phases ^ j of tiie initial 
sequence SI and the phases yQ) coming from the straight line equation Do. • 




c) The initial sequence SI is run through point by point by starting from the end (symbol of 
rank q-1); the index i is a counter that indicates the position of the phase jump (i=q-l, .-.,0). 

20 d) With each step the phases f>j (j = q - i,.. -1) are modified by + — , so that a modified 

sequence Cj^ is obtained. This sequence Ci**" is thus constituted by phases 



e) With each step one strai^t line equation Dt'*' is calculated, which straight lines determine 



where S = 2^^7. , S' 
j=o 



-2 



j.^. ,a=12/(q.(q2-l)) , ^ =6/(q.(q+l)) and 



y = 2(2q-l)/(q(q + l)) 




tbe modified sequoace C *. This equation is written as: 
25 Di^ : y=(ai)'".x + Qaif with : 
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a,+ =a.(S.')* ->9.(S, )+ and b,* =r (SX -fii^i^ 



where (S.r=^<p*=S + \.^and (S-T =Zj^j* =S'+^j(q-^i±il) 
that is to say, 

a =a +a.— a.(q = a +(A ) with (A ) = a.— j.(q ) — Bx.— 

i 0 2 2 2 o^^i^ ''i'' 2 2 2 

5 and b;=b^+ri.f -Af i.(q-^) = ^+(B;) with(B/) = yi.^-A|i.(q-i^) 

f) For each modified sequence Cj^ is calculated a mean difference {cri*)^ between the phases 
^ j of the modified sequence Ci^ and the phases y(j) coming from the straight line equation 

or. 



ij=0 

10 g) The operations c) to f)are repeated while the phases f7j (j = Q - >/" vC| - 1) of the initial 
seqiience ~ y modified. For each value of i, another modified sequence Cj is 



obtained. It is formed by die phases <p~ = 



<p_ for j = 0,...,q-i-l 



11 



^.-yfor j = ^-i,...q-l 



15 



h) Tlie final sequence S2 is formed by the sequence whose mean difference is minimal. 

For a less complex use, the mean differences are calculated based on 

the initial mean difference (c^)^ . One obtains: 



q-1 



(<r,+)2=<ro2+2«- £ -2.(B,+)5-2(Ai+)5' 
j = q-i 

+ bo [(Ai+).q.(q-l) + 2.(B|+)xi-«-j] 



(A,-^)X].(q-l).(2q-l) 



+(Bj+)x,.(q-l)-^.L(2q-i-l) 



4q.(B ; -KAi^-).(B|^).q.(q-l)f ^^■'"^^•^•^q-^^-^^^-^^ -£.(A,-^)J.(2q-,-l)-^.(B , ^)jH-^J 
' ■ 6 2 '4 

The mean difference (cr. " )^ is obtained by replacing in the expression of (cr. )^ : 



;rby-;r,(A;)by-(Ar)and(B;)by-(B,-). 
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The second example of embodiment of the means 60 for detectii^ and 
coirecting phase jun^ is described with reference to Fig. 5. In Fig. 5 is shown an example of 
an initial sequence SI comprising two phase jumps of +j . At step d) the phases are 

modified by +|- for the p symbols from rank q-k-p to q-k-1 , and by + ^ for the k symbols 

from rank q-k to q (k varies between 1 and q and p varies between 1 and q-k). 
The modified sequences obtained at step d) are thus written as: 
q>_ for j = 0,...,q-k-p-l 



15 



m +— for j = q-k-p,...,q-k-l 

j 2 

<p +7r for j = q-k,...,q 

. j 

for j = 0 q-k-p-1 

(D -— for j = q-k-p,...,q-k-l 

i 2 



q> -7C for j = q-k,...,q 
And the sttaig^it line equations calculated in step e) are written as: 
10 Dp,ic"^ : y = (ap,ic*)-x + (bp^*) with : 

(^,k)" =-(Sp^'r -A(Sp^ r and (bp^r =r.(Sp^r - ACSp^-r 
q-l 



where (Sp^ )^ = 2^ j"^ = S + p.|- + kjr and 



j-o 



j-0 

that is to say, 

a ,*=a +a[f.p.(q-k-£±i)+;ri.(q-^)]->S[pf+k;r] = a^+(A ;) 

p.k 0 2 2 z z 

with(A +) = a[f.p.(q-k-£^)+;rJc(q-^)]-;S[p^+k;r] 
P,k 2 2 2 2 

andb^/ =b^ +r[pf +k;r]-)5[|.p.(q-k-£±i)+^Jc.(q-^)] = b^ +(B^,^ 

with (B +) = y\$~+^n\- Pi^ViSi-^-Kr> + ^•k.Cq -^)] 
P.k 2 2 2 

For a less complex implementation, the mean differences (Cp^^)^ are calculated 
based on the initial mean difference {a^ . The following expression is obtained: 
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+ bo [(Apjc'').q.(q - 1) + 2Xfi^j,^).q - ^(p + 2k)] 



.q.(B^^n^ .(Ap^-).(B^^^).,.(q-l). (ApA-)^q.(q-l).(2q-l) 

- 1 -(Ap^^ )p.(2q 2k - p - 1) - ;rXAp^^ )k.{2q - k ^ 1) - .T.^^ 
The mean difference (0-^ ' )^ is derived from the expression {<^p^ ^ )^ by replacing 

Tuhy-TT, (A^/) by - (Ap fe-) and (B^ ,0 by - (B^ ,/) . 

It is evident that the means that have just been described are calculation means 
5 advantageously used in the form of a computer program intended to be executed by a 
microprocessor placed in the receiver. 

hi Fig. 6 are shown the results obtained thanks to the invention in a system as 
described with reference to Figs. 1 and 2. The curves of Fig. 6 represent the packet error rate 
(PER) plotted against the signal-to-noise ratio (SNR) in the following three cases: 
10 - curve R3 : when the phase jumps are not detected; 

- curve R2: when the phase jumps are detected and corrected with the method according to 
ttie invention; 

- curve Rl : for a Gaussian channel (no imperfection as regards either phase or frequency). 
These curves have been obtained via: 

15 - a QPSK modulation, 

- a random initial frequency error Afo comprised between -1% and +1% of the symbol 
frequency, 

- a random initial phase error 8© comprised between -^et + k, 

- packets which comprise a known preamble of 48 symbols [ai, 343] and a payload part 
20 [a49, •-.» ^592] of 544 symbols, 

- a frequency recovery algorithm (utilized by the means 23), which uses the packet 
preambles, and which is applied before the phase estimation, so that the residual frequency 
error Afi before the phase estimation is lower than or equal to 0.3% of the sjonbol frequency, 

- phase estimation means 50 which are formed by a loop which is passed through L=2 times, 
25 - on transmission, a Reed Solomon coding and a convolution coding which is punctured in 64 

states. 
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The invention is not restricted to the embodiments that have just been 
described by way of example. More particularly, it relates to any type of receiver that utilizes 
a frequency error correction device of the Tretter type and, iq)stream of this device, a phase 
estimation device based on decisions made on the received symbols. 

5 Moreover, the number of phase jumps that may be detected and corrected by 

applying this method is arbitrary. When the transmitted data packets are longer, it is possible 
to have more than two phase jumps. In that case, in order not to complicate the calculations 
too much, one advantageously chooses to divide the data packets into smaller portiona so as 
to be in the same position again where one has a maximum of two phase jumps per packet 

10 portion. The method described above is thus applied to each portion of the packet. It provides 
the position and the direction of the various phase jumps. The phases are then corrected by a 
multiple of I as a fimction of the direction of the jumps and their number. Then, the Tretter 

algorithm is applied to this corrected sequence to obtam the final frequency estimate. In order 
to avoid a jump occurring between two portions of a packet, it is desirable to provide, an 
1 5 overlap between the various portions of the same packet. 

The method proposed is generally used for any MPSK modulation by 
considering phase jumps of ± 2^ / M (instead of ± ?r / 2 for a QPSK modulation). 
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CLAIMS: 



1 . A communication system comprising at least a transmitter (1) and a receiver 

(2) intended to receive symbols coming firom a PSK modulation, and comprising estimation 
means (62) for estimating a frequency error relating to a symbol based on a sequence of 
symbol phases, 

5 characterized in that said receiver comprises calculation means (52) for calculating a phase 
sequence, called initial sequence (SI), based on decisions made on symbols, and means for 
detecting and correcting phase jxmips in this initial sequence, to supply a phase sequence, 
called final sequence (S2), to said frequency error estimation means (62). 

10 2. A communication system as claimed in claim 1 , characterized in that said 

means for detecting and correcting phase jumps comprise: 

- modifying means (100) for modifying said initial sequence (SI) so as to produce a plurality 
of modified sequences (Ci**';Cr;Cp,k"*";Cp,k^ ), which each compensate for a phase jump 
configuration, 

15 - calculation means (110) for calculating straight line equations which determine the initial 
sequence (Do) and the modified sequences ( Dj"*" ; Dj" ; Dpjc* ; Dpjc" )/ 

- calculation means (120) for calculating for the initial sequence and the modified sequences 
a mean difference between the initial or modified phases and the phases produced by the 

corresponding straight line equation ( o*^ ^ ; (o*. ^ )^ ; (cr| ~ )^ ; (cr ^ )^ ; (cXp^ " )^ ), said final 
20 sequence being formed by the sequence whose mean difference is minimal. 

3. A commimication system as claimed claim 2, characterized in tiiat said initial 

sequence is modified phase-group by phase-group. 



25 



4. A receiver intended to be used in a communication system as claimed in one 

of the claims 1 or 2. 
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5. A receiver as claimed in claim 4, characterized in that said initial sequence is 
modified phase-group by phase-group. 

6. A method of estimating a frequency error relating to a received symbol 
coming from a PSK modulation, based on a sequence of symbol phases, characterized in that 
the method comprises a calculation step of calculating a phase sequence, called initial 
sequence, based on decisions made on symbols, and a step of detecting and correcting phase 
jumps in this initial sequence, to produce a phase sequence, called final sequence, used for 
the estimation of a frequency error. 

7. A method of detecting and correcting phase jumps in an initial sequence of 
symbol phases coming from a PSK modulation, characterized in that it comprises: 

- a step (c, d) of modifying said initial sequence (SI) so as to produce a plurality of modified 
sequences which each compensate for a phase jxmp configuration, 

- a step of calculating straight line equations which determine the initial sequence (a) and the 
modified sequences (e), 

- a calculation step of calculating for the initial sequence (a) and the modified sequences (f) a 
mean difference between the initial or modified phases and the phases produced by the 
corresponding straight line equation, said final sequence being formed by the sequence 
whose mean difference is minimal. 

8. A method of detecting and correcting phase jumps as claimed in claim 7, 
characterized in that said initial sequence is modified phase-group by phase-group, 

9. A program comprising instructions for implementing the steps of a method of 
detecting and correcting phase jumps as claimed in claim 7 when said program is executed by 
a processor. 

10. A program comprising instructions for implementing the steps of a method of 
estimating a frequency error relating to a received symbol, as claimed in claim 6, when said 
program is executed by a processor. 
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